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Abstract 

Na^,KVH^ antiporters are H^-coupled cotransporters that are crucial for cellular homeostasis. Populus euphratica, a well- 
known tree halophyte, contains six NaVH^ antiporter genes {PeNHXl-6) that have been shown to function in salt tolerance. 
However, the catalytic mechanisms governing their ion transport remain largely unknown. Using the crystal structure of the 
NaVH^ antiporter from the Escherichia coli (EcNhaA) as a template, we built the three-dimensional structure of PeNHXS from 
P. euphratica. The PeNHXS model displays the typical ^M4-^M^ 1 assembly that is critical for ion binding and translocation. 
The PeNHX3 structure follows the 'positive-inside' rule and exhibits a typical physicochemical property of the transporter 
proteins. Four conserved residues, including Tyrl49, Asnl87, Aspl88, and Arg356, are indentified in the TIV14-TI\/11 1 assembly 
region of PeNHX3. Mutagenesis analysis showed that these reserved residues were essential for the function of PeNHXB: 
Asnl87 and Aspl88 (forming a ND motif) controlled ion binding and translocation, and Tyrl49 and Arg356 compensated 
helix dipoles in the TM4-TI\/11 1 assembly. PeNHX3 mediated Na^, and Li^ transport in a yeast growth assay. Domain- 
switch analysis shows that TMl 1 is crucial to Li^ transport. The novel features of PeNHX3 in ion binding and translocation 
are discussed. 
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Introduction 

Salt stress is one of the most severe environmental factors 
limiting the productivity of crop plants [1-3]. Understanding the 
mechanisms underlying plant salt tolerance, therefore, is impor- 
tant for developing crop production [4,5]. Studies have shown that 
some salt tolerant plants, named halophytes, are capable of 
surviving in saline soils [6]. Thus, exploring how the halophytes 
deal with salt is helpful to improve agricultural production [6,7]. 

Populus euphratica is a well-known tree halophyte that exists in 
saline desert areas from northwest China to western Morocco [8] . 
This tree halophyte can tolerate NaCl stress up to 450 mM [9] . Its 
seedlings accumulated low levels of Na"^ under salt stress [10]. A 
recent genomic study found that the genes encoding ion 
transporters were selectively expanded in the P. euphratica 
genome [11]. These studies suggest that maintaining low levels 
of salt in the cell is a mechanism for P. euphratica to tolerate salt 
stress. Furthermore, ion transporters controlling cellular salt 
movement may play critical roles in this tree halophyte under 
salt stress. 

Na^,K^/H antiporters (NHX antiporters) are H^-coupled 
cotransporters that transfer the Na^ or K""" across a membrane 
in exchange for protons (H''^) [12,13]. Their exchange activity is 



driven by the H''" electrochemical gradient generated by the H""" 
pumps such as the plasma membrane H'*"-ATPase or the vacuolar 
membrane H^-ATPase and H"'^-pyrophosphatase [4] . Biochemical 
and genetic studies have shown that plant NHX antiporters play 
an important role in salt tolerance [14—18]. In Arabidopsis, 
Atnhxl and Atsosl mutants are sensitive to salt stress [19,20]; 
overexpression oi AtNHXl and S0Sl/AtNHX7 reduces cytoplas- 
mic Na^ content and enhances salt tolerance in Arabidopsis 
[14,21,22]. Further study shows that SOSl activity is regulated by 
SOS2 kinase [16,17,23], and SOSl is activated by the removal of 
a C-terminal auto-inhibitory domain upon phosphorylation by the 
SOS2/SOS3 complex [24]. AtNHXl may also be regulated by 
SOS2 kinase [18]. Also, CaM binds and inhibits the Na^/H^ 
antiport activity of AtNHXl [25]. Recent study shows that 
AtNHXs may also be involved in pH and K homeostasis, vesicle 
trafficking, and plant growth and development [26—28] . AtNHX 1 
and LeNHX2 have a K"'"/H"'" transport activity and mediate K"*" 
compartmentation in vacuoles [20,29—32]. The NHX antiporters 
in Ipomoea tricolor and Ipomoea Nil are involved in vacuolar pH 
regulation; mutation of a NHX gene in Ipomea Nil abolishes the 
colour change in flowers [33,34]. nhxl nhx2 double knockout 
mutants show significantly reduced growth and abnormal stamens, 
suggesting their roles in cell expansion and flower development 
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[35]. nhxl nhx2 double mutants have reduced vacuolar K"*" pool, 
compromised turgor generation for cell expansion, and impaired 
osmoregulation, suggesting that AtNHXl and AtNHX2 are 
critical for cellular uptake and stomatal movement [36]. 

In P. euphralica, six NHX genes (PeNHXl-6) have been 
isolated and their ion transport activities have been tested in a 
yeast growth assay [37]. Using the electrophysiological technique, 
it has been found that P. euphratica has a higher Na"'"/H"'" antiport 
activity than its salt-sensitive congener P. popularis, suggesting 
that PeNHXs may function in salt tolerance in the tree halophyte 
[10]. However, understanding the catalytic mechanisms governing 
ion transport for the PeNHXs requires a structural model of these 
antiporters. Up to now, only the crystal structure of EcNhaA from 
E. Coli has been determined within the NHX gene family, and 
currently there are not any crystal structures available for the 
eukaryotic NHX antiporters in the database [38,39]. 

Homology modeling is a computational approach by which the 
three-dimensional (3D) structure of a protein (target) is constructed 
using a protein having identifi(^d crystal structure as a template 
[40,41]. Since both the bacterial EcNhaA and eukaryotic NHX 
antiporters have similar function in controlling pH and ion 
homeostasis, and share a common ancestor and a similar 
structural fold, hence, it is reasonable to use EcNhaA as a 
template to predict the structure of the eukaryotic NHX 
antiporters [42-44]. To date, EcNhaA has been used successfully 
as a template to generate the structures of Human NHXs NHE 1 
and NHA2 [43,45]. Landau et al. [43] have constructed a model 
structure of human NHEl using EcNhaA as a template, and the 
predicted structure fitted properly with the results obtained by the 
previous mutagenesis, inhibitor binding, and NMR studies [43]. 
Schushan et al. [45] generated a structure of human NHA2 based 
on EcNhaA, and some key residues involved in ion transport have 
been identified by a model guided mutagenesis analysis. The 
experimental and structural analysis revealed a novel structural 
attributes for NHA2 and suggested a mechanism of antiport 
different from the previously characterized NhaA- and NHEl -type 
transporters [45]. 

In this study, the structural features and catalytic mechanism of 
PeNHXS were studied by building a 3D structure. The structure 
of PeNHX3 was created by homology modeling. Structural 
characteristics of PeNHX3 were examined by comparing with the 
crystal structure of EcNhaA. The functions of the conserved 
residues were analyzed by mutagenesis analysis in yeast. This 
model can be used to understand the catsdytic and regulatory 
mechanisms of PeNHX3 in the tree halophyte P. euphratica 
under salt stress. 

Materials and Methods 

Fold identification, TM helix prediction, and painwise 
alignment 

The Pfam database was used to analyze protein families [46]. 
The fold- recognition server FFAS03 was used to identify protein 
folds [47]. Using the sequence of PeNHX3, FFAS03 identified the 
structure of EcNhaA (PDB ID IZCD) with a significant score of — 
121. Based on the FFAS03 analysis, scores lower than —9.5 
indicate significant similarity between the proteins; thus, EcNhaA 
is the best template for PeNHX3. 

HMMTOP [48] and TMHMM [49] were used to predict TM 
hehces. The pairwise alignments were computed with the FFAS03 
server [47] and the FUGUE server [50]. 



3D model construction by homology modeling 

The homology modeling program SwissPdb viewer was used to 
build the 3D model of PeNHXS [51]. The template EcNhaA 
(PDB ID IZCD) and the pairwise alignment beUveen PeNHX3 
and EcNhaA (Figure SI) were used to generate the structural 
model of the TM domains (residues 1 7-799). 

Yeast Strains, media, and growth conditions 

Saccharomyces cerevisiae strain AXT3 {enal-4A::HIS3, nhalA:: 
LEU2, nhxl A:: TRPl) was a gift from Dr. Jose M. Pardo [52-54]. 
This strain was the derivative of W303-1B. Untransformed strains 
were grown at 30°C in YPD medium ( 1 % yeast extract. 2 % peptone 
and 2% glucose). Transformation of yeast cells was j)crf()rmed by 
the lithium acetate method [55] . After transformation, strains were 
grown on selective HartweU's complete (SC) medium or APG 
medium (lOmM arginine, 8mM phosphoric acid, 2 mM MgS04, 
1 mM KCl, 0.2 mM CaCl2, 2% glucose, and trace minerals and 
vitamins). NaCl, LiCl, KCl, or hygromycin B was added to the 
medium. Drop test media contained 20 mM MES, and pH was 
adjusted with arginine [56] or to acidic pH values with phosphoric 
acid [57]. 

Generation of the point mutants of PeNHX3 

pDR196-PeNHX3 was used as a template to generate the point 
mutants. The site-directed mutagenesis was performed by 
Quikchange mutagenesis. The mutations were GAT to GCT, 
GAA, or AAT (D188A, D188E, or D188N); AAT to GCT, GAT, 
GAA, or CAA (N187A, N187D, N187E, or N187Q); TAT to 
GCT, GAT, AAT, GAA, or CAA (Y149A, Y149D, Y149N, 
Y149E, or Y149Q); and AGA to GCT, AAA, or GAA (R356A, 
R356K, or R356Q), respectively. For the yeast expression assay, 
these point mutants of PeNHX3 were recombined into the yeast 
vector pDR196. 

Functional expression of PeNHXS in yeast 

The AtNHX2 gene, ordered from ABRC, was cloned into the 
yeast expression vector pDR196 with the promoter PMAl. The 
empty vector pDR196 was received from Dr. John Ward as a gift. 
To clone PeNHX3, gene fragments were amplified by PCR from 
Populus euphralica cDNA using the following primers: PeNHX3 
(5'- CCGGAATTCATGGAAGCTTATATGAGTTC -3' and 
5'- CCGCTCGAGTCAATGCCATTGATTCTGAG -3'). To 
clone ScNHXl , gene fragments were amplified by PCR from the 
genomic DXA isolated from the Saccharomyces cerevisiae strain 
BJ3505 using the following primers: ScNHXl (5'-CGCGTCGA- 
CATGCTATCG4AGGTATTGC-3' and 5'-CCGCTCGAGC- 
TAGTGGTTTTGGGAAGAG-3'). The PCR fragments were 
cloned into the plasmid pDR196. AU gene fragments were verified 
by sequencing. 

AU plasmids were transformed into the yeast strain AXT3, and 
the empty vector pDR196 was transformed into the same yeast 
strains as a control. For the stress tolerance tests, yeast cells were 
normahzed in water to Agoo of 0.12. 4 |j.l ahquots of each 10-fold 
serial dilution were spotted onto AP plates supplemented with 
KCl, or YPD plates supplemented with NaCl or LiCl as indicated, 
and incubated at 30°C for 3 days. Resistance to hygromycin B was 
assayed in YPD medium. 

The salt-sensitive growth of the yeast cells expressing PeNHX3 
and its mutants was measured in liquid APG media. We prepared 
a preculture by growing the previously propagated yeast in 2 ml 
APG media overnight at 30°C. The cells were centrifuged, washed 
once with water, and resuspended in water (ODgoo~0.65). 
Inoculate 0.2 ml APG medium in a 96-weU microplate with 4 jtl 
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Figure 1. The TM segments and topology of PeNHXS. A. The TM helices were analyzed with FFAS03, FUGUE, HMMTOP, and TMHMM. The TM 
helix numbers are marked at the corresponding regions. The TM boundaries for structure modeling are highlighted in gray. B. The TM topology of 
PeNHX3. Amino acid residues are colored based on the hydrophobicity scale of Kessel and Ben-Tal (Schushan et al., 2010). 
doi:l 0.1 371/journal.pone.0104147.g001 



of the seed culture. After incubation at 30°C for 48 h, cultures 
were gently resuspended, and the ODgoo was recorded on a 
Thermo Scientific Varioskan Flash plate reader. 

Domain-switch analysis 

The AtNHXl gene, ordered from ABRC, was cloned into the 
yeast expression vector pDR196 with the promoter PMAI. A 
fragment of AtNHXl CDS beginning from the start codon ATG 
to T1092 was subcloned into the vector pDR196. Two fragments 
of PeNHXS CDS, C537 (A1102-A1638) and C399 (A1240- 
A1638), which encode two peptides of the C-terminal domain of 
PeNHXS, were fused with the AtNHXl fragment to produce the 
fusion proteins of AtNHXl-PeNHX3-C537 and AtNHXl- 
PeNHX3-C399, respectively. 

Sequence alignment analysis 

The predicted amino acid sequences of PeNHXS and AtNHXl 
were collected from the GenBank (http://www.ncbi.nlm.nih.gov). 
Proteins were compared by multiple alignments using the 
ClustalX 2.1 method (Larkin et al., 2007). The alignment is based 
on the complete amino acid sequences. 

Results 

TM helix assignment of PeNHXS 

As the only NHX antiporter having a crystal structure, EcNhaA 
has been used as a template to build the homology modelings for 
Human NHXs NHEl and NHA2 [43,45]. In this study, we are 
interested in generating a model structure for thi- P. euphratica 
NHX antiporter 3 (PeNHXS) using EcNhaA as a template. We 
found that both PeNHXS and EcNhaA belonged to the same 
phylogenetic clan as detected by the Pfam database [46] . Using the 
PeNHXS sequence as a target, we identified EcNhaA as the closest 
homologue by the fold-recognition server FFASOS [47] .These 
results suggest that both PeNHXS and EcNhaA are evolutionarily 
related and may share a similar fold. 

Since we cannot use the standard methods to align the 
sequences between EcNhaA and PeNHXS due to the low 
sequence identity between these two proteins (below 10%) (data 
not shown), we used a combined approach to obtain the alignment 
between the target and templati; sequences. The pairwise 
alignments were constructed by the FFASOS [47] and FUGUE 
algorithms [50] (Figure lA). In addition, the TM and secondary- 
structure prediction algorithms and hydrophobicity analysis were 
included to assist the selection of the TM helix boundaries 
(Figure lA). 

FFAS (the fold and function assignment system) is a profile- 
profile comparison algorithm that has a high sensitivity and 
alignment accuracy compared with other methods [47]. There- 
fore, FFAS is frequently used for the detection of remote 
homologies. In our analysis, we found that the FFASOS alignments 
match the results of FUGUE, HMMTOP, and TMHMM 
(Figure lA). Thus, we utihzed the helix assignment produced by 
the FFASOS to establish the boundaries of the TM helices for 
PeNHXS in this study (Figure lA, IB). 

Gaps were found in TM segments 2-6 of PeNHXS in the 
FFASOS alignment. We manually removed these gaps and 
assigned the helices according to the TM segments of EcNhaA 
(Figure lA). 



3D model building of PeNHXS 

To build the 3D structure of PeNGXS, we aligned the TM 
segments of PeNHXS with EcNhaA based on the FFASOS 

prediction. We then used this alignment and the EcNhaA template 
to construct a 3D model of PeNHXS with SwissPdb viewer [51]. 

As shown in Figure 2A, the PeNHXS model exhibited the 
typical TM4-TM11 assembly that has been identified in the 
structures of the NHX antiporters (EcNhaA, NHEl and NHA2) 
[38,39,43,45]. In the PeNHXS model, the TM4 and TMll 
segments unwound to form extended peptides in the center of the 
helix and crossed each other in the middle of the membrane. The 
TM4-M11 assembly was embedded between TMl, TMS, TM5, 
TMIO, and TM12 (Figure 2A). 

In the PeNHXS model, the polar residues were clustered either 
in the core structures or on extramembrane loops while the 
hydrophobic residues were facing the membrane (Figure 2A), 
displaying a typical physicochemical property of the transporter 
proteins [45]. 

The 'positive-inside' rule is an empirical observation that has 
been used to evaluate the structural model of membrane proteins 
[45]. This rule suggests that the intracellular regions are enriched 
with positively charged residues (lysine and arginine) compared 
with the extracellular regions in the structure of membrane 
proteins [45,58,59]. This 'positive-inside' rule has been well 
demonstrated by the structures of the NHX antiporters EcNhaA, 
NHEl and NHA2 [38,43,45]. We are interested in understanding 
if our PeNHXS model fits the 'positive-inside' rule. As shown in 
Figure 3, of the 25 lysines and arginines in the PeNHXS model, 18 
positively charged residues were localized in the cytoplasmic side 
whereas 7 in the extracellular side of the membrane (Figure 3), 
suggesting that our PeNHXS model follows the 'positive-inside' 
rule. 

PeNHXS mediates Na-i-, K-i- and Li-i- transport in yeast 

PeNHXS has been reported to be tolerant to salt, high K""" or 
Li""" stress in a yeast growth assay [37]. However, PeNHXS 
transport activity was not compared with the plant or yeast NHX 
antiporters. Te better understand the function of PeNHXS, we 
compared its activity with AtNHXl, AtNHX2 and ScNHXl 
using a yeast expression system. PeNHXS gene was cloned into the 
yeast vector pDR196 and expressed in a Saccharomyces cerevisiae 
strain AXTS. Strain AXTS is sensitive to salt and to high K""" since 
it lacks the functional plasma membrane Na"'"-ATPases (ScENAl- 
4), plasma membrane Na''",K''"/H''" antiporter ScNHAl, and 
vacuolar Na"'",K''"/H''" antiporter ScNHXl. The transformed yeast 
was grown on Arg phosphate (AP) or YPD medium with high 
levels of KCl, NaCl or LiCl (Figure 4). AXT3 mutants did not 
grow well in the presence of 200 mM NaCl while the nhxl- 
positive strains (W303-1B) grew dynamically (Figure 4A). The 
AXTS strain expressing PeNHX3 recovered tolerance to salt 
stress. In addition, expression of PeNHXS conferred resistance to 
high K , but to a lesser extent compared to the strains expressing 
AtNHXl, AtNHX2, or ScNHXl (Figure 4B). Interestingly, 
PeNHXS was highly tolerant to Li^ stress (Figure 4C). While 
AXT3 mutants failed to grow in the medium containing 
hygromicin B (90 |J,g/ml), PeNHXS greatly enhanced tolerance 
to the drug hygromicin B. PeNHXS tolerance to hygromicin B 
was much stronger than AtNHXl, AtNHX2, or ScNHXl 
(Figure 4D), suggesting that PeNHXS may have a significant 
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impact on endosomal compartments. These results suggest that 
PeNHX3 is similar to both the Arabidopsis and yeast NHX 
antiporters in that all NHX antiporters facilitate Na""" and K"*" 
homeostasis. However, unlike its Arabidopsis and yeast counter- 
parts, PeNHXS may play a vigorous role in Li^ detoxicity and 
endosomal trafficking. 

TM1 1 may control the Li^ transport activity of PeNHX3 

To investigate the mechanism controlling Li^ transport by 
PeNHX3, we performed the domain-switch assays. We first 
compared the sequence of PeNHXS with that of AtNHXl by 
alignment. The alignment shows that the region ranging from 
Asn368 to the C-terminal end of PeNHXS is less conserved 
compared with the remaining parts of the protein sequence (Figure 
S2). Thus, the less conserved sequence in this C-terminal domain 
suggests that the C-terminal region may be responsible for the 
diversified Li"^ transport activity between PeNHXS and AtNHXl. 
Further analysis shows that this diversified C-terminal region of 
PeNHXS mainly contains two transmembrane domains, TM 1 1 
and TM12 (Figure lA, IB). Since ion transport is primarily 
catalyzed by the transmembrane domains of the antiporters, we 
reason that these two transmembrane domains (TM 1 1 and TM 1 2) 
may be the critical players of the higher Li""" transport in PeNHXS. 
To test this hypothesis, we conducted the domain-switch analysis 
by exchanging these two transmembrane domains of PeNHXS 
with the corresponding regions of AtNHX 1 . We generated two 
constructs that produce the fusion proteins of AtNHXl and 
PeNHXS: AtNHXl-PeNHXS-C5S7 and AtNHX 1-PeNHXS- 
CS99, respectively (Figure 5A, 5B). AtNHXl-PeNHXS-C5S7 
contains a peptide ranging from Asn368 to His545 (the end of 
PeNHXS amino acid sequence), including both the TM 1 1 and 
TM12 of PeNHXS (Figxire 5B). On the other hand, AtNHXl- 
PeNHXS-CS99 contains a shorter peptide ranging from Ala414 to 
His545, including only TM12 of PeNHXS (Figure ,5B). The ion 
transport activity of these two AtNHX 1 -PeNHXS fusion proteins 
was tested in yeast (Figure 6). Interestingly, AtNHX 1-PeNHXS- 
C537 displayed a higher capacity in conferring Li""" tolerance than 
its parent protein AtNHXl, and reached the same level as 
PeNHXS (Figure 6C). In contrast, AtNHXl-PeNHXS-CS99 had 



a much lower ability in recovering Li'*' tolerance compared with 
the parent protein AtNHXl(Figure 6C). These results suggest that 
TM 1 1 may play a central role in catalyzing Li transport in 
PeNHXS. In addition, AtNHXl-PeNHXS-C5S7 conferred toler- 
ance to salt tress, high K"'" and hygromicin B, while AtNHXl - 
PeNHXS-CS99 did not (Figure 6A, 6B, 6D), suggesting that 
TMl 1 is also crucial to the transport of Na""", K""" and hygromicin B 
in PeNHXS. 

ND motif is essential for the transport activity 

Structural studies have shown that the CPA2 antiporter 
EcNhaA has two sequential Asp residues, Aspl6S and Asp 164, 
in TM5 (Figure 7A, 7B) [S8,S9]. Aspl6S and Asp 164, which form 
a DD motif, localize in the TM4-TM11 assembly region of 
EcNhaA. The DD motif has been demonstrated to be involved in 
ion binding and translocation, and thus is essential for the 
transport activity [38,S9,43,45]. In CPAl family, however, the 
corresponding residues are Asn and Asp that form a ND motif 
Study has shown that the ND motif corresponding to Asn266 and 
Asp267 in NHEl, a CPAl antiporter, has a similar function as the 
DD motif in ion binding and translocation [45]. Similar to its 
CPAl counterpart NHEl, we identified a ND motif composed of 
the residues Asnl87 and Aspl88 in the TM5 of PeNHXS, which 
are homologous to Aspl6S and Asp 164 of EcNhaA, respectively 
(Figwe 7C). To investigate the function of the ND motif in 
PeNHXS, we conducted site-directed mutagenesis by individually 
replacing Asn 187 with glutamine, glutamate, or alanine, and 
replacing Asp 188 with glutamate, asparagine, or alanine, respec- 
tively. The transport activities of the mutant genes were 
determined in a yeast expression system (Figure 8, 9). The 
transformed yeast was grown in medium with different NaCl 
concentrations or pH values (Figure 8, 9). The results showed that 
except for D188E, all the point mutants did not rescue the growth 
of the salt sensitive strain under salt stress, suggesting a specific 
requirement for an asparagine residue at Asn 187 and an acidic 
residue at Asp 188 locations, respectively. 

Studies have shown that the mutant proteins lost their transport 
activities and failed the functional complementation test when the 
DD motif in CPA2 antipoters NHA2 or EcNhaA was substituted 
to an ND motif [45,60]. We are interested in exploring if the DD 
motif is functioning in the CPAl antiporters by generating the 
N187D mutant in PeNHXS. Interestingly, the N187D mutant, in 
which the original ND motif was transformed into an DD motif, 
grew even faster than the wild type yeast under salt stress 
(Figure 8, 9), suggesting that the DD motif functions equally well 
as the ND motif in PeNHXS. 

Tyr 149 and Arg356 are vital for charge compensation 

One of the unique features of the NhaA structure is that its 
TM4 and TMll helices are discontinuous and interrupted by an 
extended chain [38, S9]. These extended chains in TM4 and 
TM 1 1 helices are flanked by partial charges formed at the N- and 
C-termini of the respective small helices connecting them. These 
partial charges are compensated by the negatively charged residue 
Asp 133 of TM4 and the positively charged residue LysSOO of 
TMIO, respectively (Figure 7A, 7B) [38,39]. In our PeNHXS 
model, we identified Arg356 in TMIO as the equivalent of LysSOO 
in EcNhaA (Figure 7C). However, we did not find a negatively 
charged residue in the position corresponding to Asp 133 of 
EcNhaA, instead we found a polar uncharged residue Tyr 149 at 
this site in PeNHXS (Figure 7C). In order to understand if these 
two residues are functioning in compensating the helix dipoles in 
the TM4-TM 1 1 assembly, we performed mutagenesis analysis 
(Figwe 8, 9). Interestingly, R356K fuUy recovered yeast growth 
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Figure 4. PeNHXS facilitate Na^ K*^ and Li^ transport and confer tolerance to Kiygromycin B in yeast. The cDNAs of PeNHXS, AtNHX2, 
AtNHX2, and ScNHXl were subcloned into the yeast expression vector pDR196 and transformed into the AXT3 mutant (enal-4 nhal nhxl). The yeast 
strains were grown overnight in AP or YPD medium. Cells were normalized in water to A600 of 0.12. Aliquos (4 |iL) from the normalized yeast cultures 
or 10-fold serial dilutions were spotted onto AP plates supplemented with KCI (B) or onto YPD plates with NaCI (A), LiCI (C), or Hygromycin B (D). The 
strains were grown at 30 C for 3 days. 
doi:1 0.1 371 /journal.pone.01 041 47.g004 



under salt stress while R356A failed to rescue salt-tolerant growth 
(Figure 8, 9), suggesting that positive charges are necessary for the 
charge compensation at the site of Arg356. Moreover, R356Q^ 
rescued yeast growth to the same level as the wild type strain 
(Figure 8, 9), suggesting that the partial positive charge of the polar 
residue glutamine can stabilize the partial negative dipoles of TM4 
and TMll in PeNHXS. 

In case of Tyr 1 49, whUe Y 1 49A did not recover yeast growth, 
all other mutants, including Y149D, Y149N, Y149E, and Y149Q, 
fully rescued yeast growth under salt stress (Figure 8, 9), suggesting 
that either the negative charges of Asparate and glutamate or the 
partial negative charges of asparagine and glutamine can stabilize 
the partial positive dipoles of TM4 and TMll in PeNHXS. 



Discussion 

PeNHXS may function in salt stress and Li-i- detoxification 
in P. euphratica 

A study has shown that PeNHXS is a tonoplast antiporter that 
mediates Na'*', K*, and Li^ transport [S7]. PeNHXS gene 
expression was highly induced in the roots, stems, and leaves of 
P. euphratica under salt stress [S7]. Our results showed that 
PeNHXS recovered yeast growth under salt stress (Figure 4A). We 
also showed that PeNHXS had a higher transport activity for Li""" 
relative to the Arabidopsis and yeast antiporters (Figure 4C), 
Thus, these results suggest that PeNHXS may play an important 
role in salt stress and Li""" detoxification in P. euphratica. 
Nevertheless, PeNHXS showed a lower K"*" transport activity 
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Figure 6. TM11 is crucial to Li^ transport in PeNHX3. The yeast 
strains were grown overnight in AP medium. Cells were normalized in 
water to A600 of 0.12 Aliquos (4 |iL) and 10-fold serial dilutions were 
spotted on AP plates supplemented with KCI (B) or YPD plates with 
NaCI (A), LiCI (C), Hyg B (D) as indicated on each panel. The strains were 
grown at 30'C for 3 days. 
doi:1 0.1 371 /journal.pone.01 041 47.g006 

compared to Na^ transport (Figure 4A, 4B), suggesting tliat 
PeNHX3 may have a higher affinity for Na""" over K"^. 

Moreover, we showed that PeNHX3 conferred resistance to 
hygTomycin B in yeast growth (Figure 4D). The change in 
sensitivity to the drug hygromycin B represents the alteration in 
membrane potential and membrane trafficking. It is increasingly 
evident that Na^,K^/H^ antiporters play an important role in 
membrane trafficking [28]. Studies have shown that yeast 
ScNlixlp is vital for vacuole fusion [61,62] and vesicle trafficking 
out of the late endosome [63-65] . Human NHE8 is critical for 
maintaining endosomal structure and for the regxilation of protein 
sorting [66]. In Arabidopsis, AtCHX17 and AtCHX21 mediate 
protein sorting [56]. AtNHX5 and AtNHX6 are involved in 
vacuolar trafficking [67]. Therefore, PeNHX3 may also play an 
important role in membrane trafficking in P. euphratica. 
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panels. The Ca atoms of the four conserved titratable residues in the 
TM4-TM1 1 assembly region are shown as spheres. A. The predicated 
structure of PeNHXS (light green) is aligned to the crystal structure of 
EcNhaA (gray) (Hunte et al., 2005). The center of the TM4-TM11 
assembly and flanking region is marked by a square. B and C. The 
marked regions of EcNhaA and PeNHX3. 
doi:l 0.1 371 /journal.pone.01 04147.g007 

The high Li^ transport activity for PeNHXS revealed in our 
study is unanticipated. We further analyzed the catalyzing 
mechanism underlying Li"^ transport in PeNHXS using the 
domain-switch technique. We obtained two fusion proteins of 
AtNHXl and PeNHXS by exchanging two peptides from the C- 
terminal domain of PeNHXS with the corresponding regions of 
AtNHXl (Figure 5 A, 5B). The results showed that the fusion 
protein AtNHXl-PeNHX3-C537, which contains the TMl 1 and 
TM12 of PeNHXS, significantly improved yeast growth under Li"*" 
stress, whereas the fusion protein AtNHXl-PeNHX3-C399 that 
contains only TM 1 2 did not show any capability in conferring Li 
tolerance in yeast (Figure 6C). Moreover, we found that AtNHXl - 
PeNHXS-C5S7 also improved tolerance to salt stress, high K"^ and 
hygromicin B, but AtNHXl-PeNHX3-C399 did not (Figure 6A, 
6B, 6D). Taken together, these results suggest that the TMl 1 may 
play an important role in catalyzing Li^ and other ion transport in 
PcXIIXS. However, our results are preliminary, and more studies, 
including site mutational analysis, are needed to understand the 
mechanism controlling Li transport in PeNFIXS in yeast. In 
addition, it will be interesting to investigate the in vivo function of 
PeNHXS in Li transport and detoxification by complementing 
Arabidopsis mutants under stress conditions. 



The TM4-TM11 assembly is critical for maintaining the 
structure and ion transport in PeNHXS 

The unique feature of the crystal structure of EcNhaA is that 
there is an assembly of the TM4 and TMll segments lying in its 
structural core, where the two helices cross each other in the 
middle of the membrane [S8]. In the TM4-TM11 assembly, the 
helices are unwinding to form extended peptides in the center of 
the helix, providing structural basis for the ion-binding site [38]. 
The assembly forms dipoles that are charge-compensated and 
stabilized by Asp 133 in TM4 and LysSOO in TMIO, respectively 
[38] (Figure 7A, 7B). In addition, the DD motif containing Aspl63 
and Aspl64 in TM5, which localizes in the TM4-TM1 1 assembly 
region of EcNhaA, is involved in ion binding, translocation, and 
transport activity (Figure 7A, 7B) [38,39]. Interestingly, both the 
model structures of NHEl and NHA2 contain a similar TM4- 
TM 1 1 assembly as well as the conserved residues in their 
structural cores [43,45]. Mutagenesis analysis demonstrated that 
the conserved residues in both NHEl and NHA2 were vital for ion 
binding and partial charge compensation [43,45] . In our PeNHXS 
model, we also indentified a TM4-TM1 1 assembly that contained 
extended peptides in their helices (Figure 7C). Four conserved 
residues, Asnl87 (in TM5), Asp 188 (in TM5), Tyr 149 (in TM4) 
and ArgS56 (in TMIO), which are the equivalents of Asp 163, 
Asp 164, Asp 133, and LysSOO in EcNhaA, respectively, have been 
identified in the PeNHXS model (Figure 7). Mutagenesis analysis 
showed that these four conserved residues were essential for ion 
transport in PeNHXS (Figure 8, 9). Thus, these results suggest that 
the TM4-TM1 1 assembly, a structural core common to the NHX 
antiporters, may play important roles in ion binding and transport 
in PeNHXS. 

An alternating access mechanism of ion transport has been 
suggested for EcNhaA [38,39]. In this hypothesis, EcNhaA forms 




Figure 8. Functional analysis of the conserved residues in the TiVI4-TM1 1 assembly region of PeNHX3 (dosage response). The cDNAs 
of the wild type PeNHXS and the PeNHXS mutants were subcloned into the yeast expression vector pDR196 and transformed into the AXTS mutant 
(ena1-4 nhal nhxl). Yeast cells were grown in APG media containing different concentrations of NaCI. Yeast growth was determined at 600 nm 
(ODsoo) after 48 h at SO°C. Data are the average of three experiments. Asterisks indicate significant difference (PsO.05; f test). 
doi:1 0.1 371/journal.pone.0104147.g008 
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two funnels that lay the path for ion-transport. One funnel opens 
to the cytoplasm while the other to the periplasm. However, these 
two funnels do not form a continuous pore, but are separated by 
densely packed non-polar residues in the middle of the membrane. 
Thus, conformational alterations are required to switch the ion- 
binding site betwc'(;n inward and outward orientations. The TM4- 
TMIl assembly may provide the structural basis for the 
alternating access mechanism of ion transport, in which the 
alternating conformation of the ion binding site is induced at the 
bottom of the funnels upon pH activation [38,39]. Since PeNHX3 
displays a similar TM4-TM1 1 assembly in its model structure, we 
thus suggest that PeNHXS may catalyze its ion transport with the 
alternating access mechanism. 

Asn187 and Asp 188 forms a ND motif that is crucial for 
ion binding and ion translocation in PeNHX 

Studies have shown that the DD motif of CPA2 NHXs or the 
ND motif of CPAl NHXs plays an important role in ion transport 
[43,45]. The DD or ND motif, which is localized close to the 
TM4-TM 1 1 assembly of the NHX structures, is involved in cation 
binding and thus is essential to the ion transport of NHX 
antiporters [43,45]. Our mutagenesis studies showed that the ND 
motif of PeNHXS, which is composed of the residues Asnl87 and 
Asp 188 in the TM5, was essential for yeast grow th under salt stress 
(Figure 8, 9), suggesting that the ND motif of PeNHXS plays an 
important role in ion transport. 

An interesting question to ask: is the DD or ND motif 
exchangeable in CPAl and CPA2 antiporters? Studies have 
shown that in NHA2, a CPA2 NHX, the ion transport activity was 
completely inhibited when its DD motif was mutated into ND 
motif [45], suggesting that the DD motif is strictly required for 



CPA2 NHXs. However, our experiments showed that yeast 
growth was not affected when the ND motif of PeNHXS was 
replaced by a DD motif (Figure 8, 9), suggesting that the DD motif 
is fully functional in CPAl NHXs. 

Since both the DD and ND motifs work equally well in CPAl 
NHXs, then why CPAl NHXs retain a ND motif but not a DD 
motif? Or in other words, what is the function of the ND motif in 
CPAl NHXs? Landau et al. [43] proposed that the ND motif in 
CPAl controls the Na''":H"'' stoichiometry of their transport 
activities. They suggested that Aspl64 in EcNhaA or Asp267 in 
NHEl serves to alternately bind Na^ or H"^. While Asp 163 in 
EcNhaA binds another proton, its corresponding residue Asn266 
in NHEl is not involved in cation binding. Asn266 may function 
in maintaining the structure of the antiporter. Therefore, 
mutagenesis studies are needed to indentify the function of the 
Asp and Asn residues in controlling the Na"'":H"'" stoichiometry in 
CPAl or CPA2 antiporters. 

PeNHXB has a novel feature for the compensation of 
helix dipoles in the TM4-TM11 assembly 

In NhaA, the extended chains in TM4 and TM 1 1 helices are 
flanked by partial charges [38,39]. These partial charges are 
stabilized by the negatively charged residue AsplSS and the 
positively charged residue LysSOO (Figure 7A, 7B) [38,39]. In 
PeNHXS, we identified residues Tyrl49 and Arg356 that were 
corresponding to AsplSS and LysSOO of EcNhaA, respectively 
(Figure 7C). Our mutagenesis analysis showed that Tyrl49 and 
ArgS56 were also fimctioning in compensating the partial charges 
of the helix dipoles in the TM4-TM1 1 assembly in PeNHXS 
(Figure 8, 9). 
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Interestingly, a negatively charged residue that is equivalent to 
LysSOO in EcNhaA exists in all NHX antiporters with a model 
structure, including NHEl (Arg425), NHA2 (Arg432) and 
PeNHX3(Arg356) (Figure 7; [43,45]), suggesting that this nega- 
tively charged residue is highly conserved within NHX gene family 
and play essential roles in ion transport. On the other hand, the 
negatively charged residue corresponding to Asp 133 in EcNhaA is 
less conserved (Figure 7; [43,45]). So far, Asp238 in NHEl is the 
only residue that has been identified to be the equivalent of 
Asp 133 of EcNhaA [43]. PeNHX3 has a polar non-charged 
Tyrl49 (Figure 7C), and NHA2 does not have a negatively 
charged residues corresponding to Asp 133 [45], suggesting that 
there are variations at this site for charge compensation. This 
implies that PeNHX3 and NHA2 may have a novel mechanism 
for charge compensation. Indeed, in NHA2, the partial positive 
charge of the helix dipoles can be stabilized by the hydroxyl side 
chains of Ser245 and Thr462 at the N termini of the helices 
connecting the ("xtended chains of TM4 and TM 1 1 , respectively 
[45]. Thus, more studies are needed to identify the restudies 
involved in charge compensation in PeNHXS. 

In conclusion, the PeNHX3 structure created by homology 
modeling in this study sheds light on its function and catalytic 
mechanism. Our results demonstrate that PeNHX3 has unique 
features in ion binding and translocation while sharing the same 
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Structure fold and catalytic mechanism as its bacterial and human 
counterparts. This model structure may facilitate understanding 
the function and regulatory mechanism of PeNHX3 in the tree 
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